Instability of diverging and converging flows in an annulus 
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Abstract 

The stability of two-dimensional diverging and converging flows in an annulus between two 
permeable cylinders is examined. The basic flow is irrotational and has both the radial and 
azimuthal components. It is shown that for a wide range of the parameters of the problem, 
the basic flow is unstable to small two-dimensional perturbations. The instability is inviscid 
and oscillatory and persists if the viscosity of the fluid is taken into consideration. 



1 Introduction 



Or 



In this paper we present a simple example of an instability of a steady inviscid flow in an 
annulus between two permeable circular cylinders. Despite the utter simplicity of the basic flow, 
this instability appears to be overlooked by other researchers. The basic flow is irrotational 
and rotationally symmetric. Both the radial and the azimuthal components of the velocity are 
inversely proportional to the radius. We show that if the ratio of the azimuthal component of the 
velocity to its radial component is larger than a certain critical value, then the flow is unstable to 
small two-dimensional perturbations. The instability is oscillatory: the neutral modes represent 
azimuthal travelling waves. The most striking feature of the instability is that it is present 
not only in the diverging flow but also in the converging flow. This contradicts the common 
view that diverging flows are always unstable while converging f lows are stable, which originates 



in the studies of the Jeffer y-Hamel flow (e.g., Goldshtik et al . 1991 ; Shtern Hussain . 19931 ; 



McAlpine fc Drazin I . ll998T ) 



The stability of viscous flows between permeable rotating cylin ders with a radial flow to three- 
dimensional perturbations had been studied by many authors (e. g., Bahl L 197C ]_ Chang &i Sartorvl . 



19671 : iMin fc Lueptow 1 . 11994 ; iKolvshkin fc Vaillancourt I . Il997l ; iKolesov & Shapakidze I . Il999l ). 



One of the main aims of these studies was to determine the effect of the radial flow on the 
stability of the circular Couette- Taylor flow, and the general conclusion was that the radial flow 
changes the stability properties of the flow: both a converging radial flow and a sufficiently strong 
diverging flow have a stabilizing effect on the Taylor vortex instability, but when a divergent 



flow is weak, it has a destabilizing effect (|Min fc Lueptow I . 11994 : IKolvshkin &; Vaillancourt 
19971 ). Since the above papers had been focused on the effect of the radial flow on the known 



instability, it remained unclear whether the radial flow itself can induce instability for flows 
which are stable without it. The results described below may help to answer this question. It 
should be mentioned that the purely azimuthal flow with the veloci ty inversely proportion al to 



the radius is stable not only to two-dimensional perturbations (see iDrazin &: Reid I . Il98ll ) but 



also to three-dime nsional perturbations (th is can be deduced from the sufficient condition for 
stability given by iHoward fc Gupta I (jl967T )h The results of our study show that this stable 
flow becomes unstable if the radial flow is present and the weaker this radial flow is, the more 
unstable the flow becomes. 

The outline of the paper is as follows. In Section 2, we formulate the problem. Section 
3 contains the inviscid stability analysis of the diverging flow. In Section 4 we show that the 
inviscid stability results of Section 3 represent a valid asymptotic of the corresponding viscous 
problem in the limit of high Reynolds numbers. In Section 5, we briefly discuss the stability of 
the converging flow. Finally, conclusions are presented in Section 6. 
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2 Formulation of the problem 



We consider two-dimensional inviscid incompressible flows in an annulus D between two concen- 
tric circles with radii r\ and r<i {r% > T\). The circles are permeable for the fluid and there is a 
constant volume flux 2ttQ of the fluid through the annulus (the fluid is pumped into the annulus 
at the inner circle and taken out at the outer circle). Suppose that r% is taken as a length scale, 
r\jQ as a time scale, Q/r\ as a scale for the velocity and pQ 2 /rf for the pressure where p is the 
fluid density. Then the two-dimensional Euler equations, written in non-dimensional variables, 
have the form 
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Here (r, 9) are the polar coordinates, u and v are the radial and azimuthal components of the 
velocity and p is the pressure. It is known that if there is a non-zero flow of the fluid through 
the boundary, it is necessary to prescribe additional boundary conditions on the part of the 
boundary where the fluid enters the flow domain. What conditions should be added is a subtle 
question and there are s everal answers that lead to mathematically correct initial boundary 



value problems (see, e.g., lAntontsev et al I . Il990l : iMorgulis &: Yudovich I . l2002f ). We will use the 



boundary condition for the tangent component of the velocity, which at the first approximation 



corresponds to the condition at a porous cylinder (see iBeavers &: Joseph 1. 119671) and for which 



the corresponding mathematical problem is well-posed (e.g., Antontsev et al . 1990l ). So, our 
boundary conditions are 



1, 



1/a, 



=1=7, 



(4) 



where a = Tijr\ and 7 is a constant that represents the ratio of the azimuthal velocity to the 
radial velocity at the inner circle. 

Problem ([TJ-Q has the following simple rotationally-symmetric solution: 



u(r, 9) = 1/r, v(r, 9) = 7/r. 
In the next section we investigate the stability of this steady flow. 



(5) 



3 Inviscid stability analysis 

We consider a small perturbation (u,v,p) in the form of the normal mode 

{u,v,p} = Re \{u{r),v(r),p(r)}e at+me 
where neZ, This leads to the eigenvalue problem: 
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and 

u(l) = 0, u(a) = 0, v(l) = 0. (10) 

First we note that the eigenvalue problem ([7|)- (|10p has no nontrivial solution for n = 0. Indeed, 
Eq. U for n = and the boundary conditions for u imply that u = 0. Equation (jHj) yields 
v = C r~~ l e~ ar I 2 where C is an arbitrary constant. Substitution of this into the boundary 
condition v(l) = leads to the conclusion that C = 0. 

Now consider eigenvalue problem ((7|)- (jl0j) for It is convenient to introduce the stream 

function tp(r) such that 

u = — 4>(r), v = —ip'(r). 
After elimination of the pressure, we obtain 

$(l) = ^(a) = 0, V'(1)=0, (12) 
1 - n 2 . 

= *l>" + - (13) 

L$ = Ce~ h{r) (14) 
where C is a constant and h(r) = ar 2 /2 + in7lnr. The general solution of (|14|) is 

r 

^ = — + C 2 r n + — / (r n s- n+1 - r- n s n+1 ) e~ h{s) ds (15) 
r n 2n J 



where 



It follows from (fTTl) that 



where C\ and C2 are arbitrary constants. Substitution of (jlojl into the boundary conditions (|12j) 
yields a system of linear equations for Ci, C2 and C. The requirement that this system has a 
non-trivial solution results in the dispersion relation for a: 

a 

D(a, n, a, 7) = / (a n r - n+1 - a - n r n+1 ) e — 2 /2-m 7 inr dr = Q (lg) 



A few conclusions can be made just by looking at the dispersion relation. First, for the 
purely radial flow (7 = 0), it shows that (i) there are no real eigenvalues (because the integrand 
in (fTB"j) is non-negative for all r £ [l,o]) and (ii) the eigenvalues appear in complex conjugate 
pairs (i.e. if a is an eigenvalue, so is its complex conjugate a). Second, it is easy to see that 



D(a,n,a,-y) = -D(cr, -71,0,7), (17) 
D(a, n, a, 7) = —D(a, — n, a, —7). (18) 

These relations imply it suffices to consider only positive n and 7. 

Numerical evaluation of (|16j) shows that (i) the purely radial flow (7 = 0) is stable for 
any a > 1, i.e. Re(er) < for all azimuthal modes, and (ii) for each azimuthal mode, when 
7 increases, the egenvalues in the upper half plane (Im(cr) > 0) move to the right and the 
egenvalues in the lower half plane (Im(<r) < 0) move to the left, and there is a critical value 
7 C > of parameter 7 at which one of the eigenvalues crosses the imaginary axis, so that 

Re(er) > for 7 > 7 C and Re(cr) < for 7 < j c . 
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Figure 1: Neutral curves for n = 1, . . . , 6. The region above each curve is where the correspond- 
ing mode is unstable, (a) - diverging flow, (b) - converging flow. 
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Figure 2: Azimuthal wave number n of the most unstable mode as a function of a. (a) - diverging 
flow, (b) - converging flow. 



Neutral curves (Re(cr) = 0) on the (a, 7) plane are shown in Fig. QJa) for n = 1, . . . , 6. For each 
azimuthal mode, the critical circulation j c is a decreasing function of a and tends to a limit 
value as a — > 00. The calculations also show that, for each n, 7 C — >■ 00 as a — >■ 1. For large a, 
the first mode that becomes unstable when 7 increases from (we will call it the most unstable 
mode) has the azimuthal wave number n = 2. When a decreases, first the mode with n = 3 
becomes the most unstable one, then the mode with n = 4, and so on. The wave number n 
of the most unstable mode versus a is shown in Fig. [2^ a). The jumps in n correspond to the 
intersection points of the neutral curves in Fig. [TJ 

This instability is oscillatory. The neutral modes represent azimuthal waves travelling with 
the phase speed A = — Im(<j) /n in the counterclockwise direction. The graphs of A versus a for 
n = 1, . . . , 6 are shown in Fig. [3]^ a). For n = 2, . . . , 6, the phase speed is positive for all a, while 
for the neutral mode with n = 1, it changes sign: it is positive for relatively small a and negative 
(though small in magnitude) for large a. 
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Figure 3: The phase speed A = —Im{a)/n of the neutral modes versus a for n = 1, . . . , 6. (a) 
diverging flow, (b) - converging flow. 



4 Effect of viscosity 

The steady flow @ is also the solution of the Navier-Stokes equations that satisfies (HJ) and the 
additional boundary condition 

v\ r=a =l/a. (19) 

This solution corresponds to the situation when the inner and outer cylinders rotate anticlockwise 
with angular velocities 7 and 7/a 2 respectively. The aim of this section is to show that for 
sufficiently high Reynolds numbers the unstable inviscid modes found in section 3 give a good 
approximation to the corresponding viscous modes. 
The viscous counterpart of Eq. (fTT|) is 



a + —± + - d r 



Li> = ±- L 2 ip, 
Re 



(20) 



where Re = Qjv is the Reynolds number (y is the kinematic viscosity of the fluid). Equation 
(f2"0"|) and the boundary conditions (cf. (fT2"j) ) 



■0(1) = 0, ^'(1) = 0, ^(o) = 0, ip'(a) = 



(21) 



represent an eigenvalue problem for a. In what follows we will construct an asymptotic ex- 
pansion of the solution to this eigenvalue problem in the limit Re — > 00. It is known (e.g., 



. 2000; 


Yudovich . 


2001: 


Ilin . 



lem like this involves a boundary layer at the part of the boundary where the fluid leaves the 
domain. So, we assume that the asymptotic expansion has the form 



a = o"o + Re Vi + Re 2 o\ + . . . , 

V> = ^ (r) + fle -1 $i(r) + Mv)] + Re- 2 [^ 2 {r) + Mv)} 



(22) 
(23) 



Here rj = (a — r)Re in the boundary layer variable. Functions ipk{ r ) (k = 0, 1, . . .) represent the 
regular part of the expansion, and 4>k(v) (^ = 0, 1, . . .) give us boundary layer corrections to 
the regular part. The latter ones are needed to satisfy the additional boundary condition that 
appears in the viscous problem. We assume that the boundary layer part rapidly decays outside 
thing boundary layer near r 



a: 



as rj — y 00, 



(24) 
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for every s > and for each k = 0, 1, . . . 

The regular part of the expansion is obtained by substituting {22]) and the formula 

i> = ip (r) +e$i(r) + e 2 4(0 + . . . (25) 

into Eq. (j20|) . This yields 

+ ^ + - ^ ^4 = 0, (26) 

cr H =- + - <9 r L^i = —a\LipQ + LLi/j , etc. (27) 

r z r / 



To obtain boundary conditions at r = 1, we substitute (I25p into the first two boundary conditions 
PUD . This yields 

4(1) = ^(1) = (28) 

for each = 0, 1, . . . For the boundary conditions at r = a, the boundary layer part must be 
taken into account. So we substitute (|23p into the last two conditions (|21 j) and use the relation 
d y = —Re d v . As a result, we get 

4(a) = 0, ^(a) - 0' O (O) = 0, (29) 
4(a) + 00 (0) = 0, 4(a) - 0'i (0) = 0, etc. (30) 

Comparing Eqs. {25]), {S3]) (for fc = 0) and the first of the conditions ([29]) with {UJ and {SJ, we 
conclude that, at the leading order, we have the inviscid eigenvalue problem. Note that at this 
stage we cannot satisfy the second of the conditions {29]), and that is why we need a boundary 
layer at r = a. 

To derive equations for the boundary layer part of the expansion, we substitute {23]) into Eq. 
(|20|) and take into account that V'fc (k = 0,1, . . .) satisfy (|26[) . {27]). Then we make the change 
of variables r = a — Re^ 1 r/, expand every function of a — Re -1 rj in Taylor's series at Re~ 1 = 
and, finally, collect terms of the equal powers in e. In the leading order, we obtain 

-i^0o = ^0o. (31) 

This should be solved subject to the condition of decay at infinity and the second condition {29]) 
that can be written as 

O (O) = 4(a). (32) 

The solution of (I3TD. ([32]) is 

0o (r?) = -a4(a) e -^. (33) 

Note that the boundary layer does not affect the leading order eigenvalue do- It can be shown 
that the boundary layer will affect the first order viscous correction a\ , but we will not compute 
it here. 



5 Converging flow 

In the previous two sections we have found that the diverging flow {5]) can be unstable in the 
framework of the inviscid theory and that the instability persists if the viscosity is taken into 
consideration. The natural question to ask is whether this instability occurs only in diverging 
flows or a similar converging flow is also unstable? To answer this question we consider the 
converging flow 

u{r,6) = —l/r, v(r,9)=j/r, (34) 



6 



which differs from ([5]) only by the sign of the radial component of the velocity. This is a 
rotationally-symmetric solution of Eqs. CE])-© that satisfies the boundary conditions 

u\ r=1 =l, u\ r=a =^, v\ r=a =~f. (35) 

Note that now the third boundary condition (for the azimuthal component of the velocity) is 
imposed at r = a (i.e. on the part of the boundary where the fluid flows into the domain). 

A repetition of the analysis of Section 3 leads to the following dispersion relation for the 
eigenvalue a: 

a 

D!((T,n,0,7) = J (r~ n+1 - r n+l ) e ^ 2 /2+inylnr dr = Q _ ( 3 g) 

1 

This dispersion relation is very similar to the one we had before (cf. (|16h ) and satisfies the same 
symmetry relations (|17p and (|18p . The neutral curves in the a -7 plane are shown in Fig. QJb) 
for the modes with n = 1, . . . , 6. The wave number n of the most unstable mode as a function 
of a is shown in Fig. EJb). Unlike the diverging flow, here the mode with n = 1 is the most 
unstable one for sufficiently large a. The phase speed A = — Im(cr)/n is shown in Fig. [3]^b). In 
contrast with the diverging flow, the phase speed for the converging flow is positive for all n > 1 
and its magnitude is almost the same for all modes with n > 3. 



6 Conclusions 

We have shown that a simple inviscid irrotational flow between two permeable cylinders are 
unstable to small two-dimensional perturbations. The instability is oscillatory and persist if the 
viscosity of the fluid is taken into account. It is a genuinely new instability in the sense that 
there is no instability if either the radial flow or the azimuthal flow is absent. 

As we discussed in Section 4, the inviscid flow ([5]) is also an exact solution of the Navier- 
Stokes equations corresponding to a particular case of the flow between permeable rotating 
cylinders when = r\jr\. In the general steady rotationally symmetric viscous flow, 

the radial component of the velocity is the same as before, while the azimuthal component is 
v = Ar 1+ ^ e + Br" 1 where A and B are constants depending on 71 = VL\r\/Q, 72 = ^r^/Q, 
a = 7^2 /n and Re. In the limit of high Reynolds numbers this flow is well approximated by ([5]) 
everywhere except for a boundary layer at the flow outlet. This fact suggests that the instability 
described above may be relevant to the whole family of rotationally symmetric viscous flows. 
This, nowever, requires a further analysis and is a topic of a continuing investigation. 

There are many interesting question that are left unanswered in the present paper. For 
instance, it is known that in problems with nonzero flow through the boundary, the stabil- 
ity may depend on the bo undary conditions imposed at the inlet and the outlet (see, e.g., 
Gallairea fc Chomaz The question that arises in this context: will the same flow be 



unstable if different boundary conditions are imposed? Another interesting question: what are 
the stability properties of the flow ([5]) with respect to three-dimensional perturbations? At the 
moment, these are open problems for future investigations. 

We are grateful to Professor V. A. Vladimirov for helpful discussions. 
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